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I. INTRODUCTION 


In the last two decades, the NRQCD color octet mechanism [1] for heavy quarkonium 
production has been very successful in explaining experimental data at high energy colliders 
such as at Tevatron j^J and at LHC jsj]. In its original formulation [lj the proof of factoriza¬ 
tion in heavy quarkonium production in NRQCD color octet mechanism was lacking. The 
proof of factorization is an essential requirement to study heavy quarkonium production 
at high energy colliders. Factorization refers to separation of short-distance effects from 
long-distance effects in quantum held theory. 

Recently the proof of factorization in heavy quarkonium production in NRQCD color 
octet mechanism is given at next-to-next to leading order (NLLO) in coupling constant by 
using diagrammatic method of QCD 4-6]. However, the proof of factorization in heavy 
quarkonium production in NRQCD color octet mechanism at all order in coupling constant 
is still missing. In this paper we will prove factorization in heavy quarkonium production 
in NRQCD color octet mechanism at all order in coupling constant by using path integral 
method of QCD. 

The typical non-perturbative NRQCD matrix element in heavy quarkonium production 
is given by [11] 


< o|ego >=< o| x t (o)^(o)(aW)e t (o)R:;x(o)|o > 


(i) 


where £ is the two component Dirac spinor held that annihilates a heavy quark, y is the two 
component Dirac spinor held that creates a heavy quark, a) H is the operator that creates the 
heavy quarkonium H in the out state. The factors K n and K' n are products of a color matrix 
(either a unit matrix or T a ), a spin matrix (either a unit matrix or a 1 ) and a polynomial of 
covariant derivative D. The color and spin indices on the helds y and £ have been suppressed. 

The production cross section for heavy quarkonium H at transverse momentum Pt in 
NRQCD factorizes into a sum of perturbative functions times universal matrix elements, 

d<7A+B^H+X(P T ) — Y ddA+B^QQM+X(Prr) <^n> ( 2 ) 

n 

where each NRQCD non-perturbative matrix element < O n > represents the probability of 
a heavy quark-antiquark pair in state [n], such as color singlet or color octet etc., to produce 
heavy quarkonium state H. 
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The fragmentation function for parton i to evolve into a heavy quarkonium at large Pt 
is factorized according to [7] 


D H /i(z, m c , n) = d i^QQ[n](z, rn c , /i) < O n > 


( 3 ) 


in terms of same NRQCD non-perturbative matrix elements, along with perturbative func¬ 
tions d^QQ^z, m c , /i) that describe the evolution of an off-shell parton into a heavy quark- 
antiquark pair in state [n], such as color singlet or color octet etc.. 

At a first glance it can be easily seen that the non-perturbative NRQCD matrix element 
in eq. (JI]) is not gauge invariant unless it is a color singlet S-wave non-perturbative matrix 
element. Hence one expects that any non-canceling infrared divergences in the perturbative 
Feynman diagrams of heavy quark-antiquark production short-distance coefficient can not 
be factorized in to the definition of the non-perturbative NRQCD matrix element in eq. © 
to study heavy quarkonium production at high energy colliders in the NRQCD color octet 
mechanism. 

This is explicitly shown in jJ-6] where the NNLO coupling constant calculation shows 
that the above non-perturbative NRQCD matrix element in eq. © is not consistent with 
factorization of infrared divergences unless it is a color singlet S-wave non-perturbative 
matrix element. By using the calculation at NNLO in coupling constant and to all powers in 
the heavy quark relative velocity it was shown in [d-6| that the octet S-wave non-perturbative 
NRQCD matrix element which is gauge invariant and is consistent with the factorization of 
infrared divergences is given by 


< 


0|a|0 >=< 0| X t (0)AA^(0)# )t (0) e6 (al,aH)# ) (0) fea e t (0)<, a x(0)|0 


> 


( 4 ) 


where 

$! A) (0) = VeM~igT {A)c dXl • ( T {A)c ) ab = -*/ abc (5) 

Jo 

is the gauge link or the non-abelian phase in the adjoint representation of SU(3), A^ a {x) is 
the gluon field, V is the path ordering and P is the light-like four-velocity. 

Note that a necessary condition for NRQCD factorization is that the long-distance behav¬ 
ior of the non-perturbative NRQCD matrix element must be independent of the light-like 
vector P. Such a dependence would be inconsistent with NRQCD factorization because 
the infrared divergences of < O n > must match those of cross sections, in which there is 
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no information on ZC In ^-6| we have verified the Z M independence of the infrared pole at 
NNLO in coupling constant and to all powers in heavy quark relative velocity. 

Since the NRQCD matrix element in eq. (IT]) is a non-perturbative quantity in QCD it 
can not be calculated by using perturbative QCD methods. It is well known that a non- 
perturbative function can not be studied by using perturbative methods, no matter how 
many orders of perturbation theory is used. Hence path integral formulation (as opposed to 
diagrammatic method using perturbation theory) is necessary to study properties of non- 
perturbative quantities in QCD at all order in coupling constant. The only path integral 
formulation to study factorization of soft and collincar divergences at all order in coupling 
constant in quantum field theory is given by R. Tucci in [8] which is exact for QED but 
is not exact for QCD. We have extended this path integral approach to QCD to prove 
factorization in QCD at all order in coupling constant in |9(. In this paper we will extend 
this path integral approach to prove NRQCD factorization at all order in coupling constant 
in heavy quarkonium production. We will prove that the long-distance behavior of the non- 
perturbative NRQCD matrix element is independent of the light-like vector Z M at all order 
in coupling constant. 

The paper is organized as follows. In section II we briefly describe the lagrangian density 
in NRQCD and in QCD. In section III we discuss infrared divergences in NRQCD and 
in QCD. In section IV we include heavy quark in the path integral formulation of QCD. 
In section V we describe infrared divergences in NRQCD and the light-like Wilson line in 
QCD. In section VI we show that the eikonal current of the light-like charge generates pure 
gauge field in quantum field theory. In section VII we show how the pure gauge field in 
quantum field theory can be used to describe soft (infrared) divergences. In section VIII we 
study heavy quark-antiquark non-perturbative matrix element in the presence of light-like 
Wilson line in QCD. In section IX we prove factorization in heavy quarkonium production 
in NRQCD color octet mechanism at all order in coupling constant and to all powers in the 
heavy quark relative velocity. In section X we show that the factorization theorem is a key 
ingredient in calculation of NRQCD heavy quarkonium production cross section at collider 
experiments. Section XI contains conclusions. 
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II. LAGRANGIAN DENSITY IN NRQCD AND IN QCD 


The lagrangian density in QCD including heavy quarks is given by 10] 


£qcd = -- A F° lu F> M ' a + £ U'fD,,. - mMn + - M]T 

4 z=i 


( 6 ) 


where F^ ua is the full non-abelian gluon field tensor, ipi is the Dirac field of the light quark 
(/ = u,d,s), T is the Dirac field of the heavy quark, ZT‘ is the covariant derivative, is the 
Dirac matrix, mi is the mass of the light quark and M is the mass of the heavy quark. 

In NRQCD an ultraviolet cutoff A ~ M is introduced. The lagrangian density in NRQCD 
is given by jlj 


£nrqcd — 




light 


C 


heavy 


6 C 


(7) 


where 

flight = —F^F*™ + y i>,[D - mi]ipi, ( 8 ) 

4 1=1 

•7-heavy = f f [»A + + xH iD t ~ 7^j]X (9) 

and 

^ = ^ [ e t(D2)2 e-x t(D2)2 A ] 

+ ^[e t(D -gE-gE- D)£ - X t( D -gE-gE- D) X ] 

xjE-jEx iD) • at - x f (xD xjE-jEx xD) • a X ] 
+^[^gE-crt-X^gE-ax] +■■■ ( 10 ) 

where D t and D are the time and space components of the covariant derivative D M and E 
and B are electric and magnetic components of the gluon field tensor and cr is the Pauli 
spin matrix. The dimensionless coefficients ci, C 2 , C 3 , C 4 etc. in eq. (ITU are obtained by 
matching NRQCD with QCD ]lj. 


III. INFRARED BEHAVIOR IN NRQCD AND IN QCD 


Note that in order for the factorization formula to hold in eqs. (J2]) and ([3]) the perturbative 
functions have to be infrared-safe by definition because infrared limit corresponds to long¬ 


distance regime jl]. However, as found 
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5] the NNLO infrared pole contribution to 
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order v 2 is given by 


r?2 


E (P, v, l ) = a\ 


1 v 
3eT 


( 11 ) 


which is not zero where v is the relative velocity of the heavy quark-antiquark pair. Eq. CD 
is in the rest frame of the heavy quarkonium (P = 0) where P M is the four-momentum of the 


heavy quarkniuom and P is the four-velocity of the 


ight-like Wilson line which is fixed to be 
6 ]. The presence of non-zero infrared pole 


P = <P _ along the minus light cone direction in [ 4 J- 6 I 
in eq. (fill) implies that infrared poles will appear in perturbative functions at NNLO and 
beyond when the factorization is carried out with octet non-perturbative NRQCD matrix 
element < ’x}K n ti(a' H aH)^K' n x > in the conventional manner as given by eq. (HD in eqs. (J2D 
and ©. On the other hand, when defined according to its gauge-completed form as given 
by eq. (J4J) each octet non-perturbative NRQCD matrix element itself generates precisely 
the same pole terms given in eq. dill) above. This conclusion is valid to all powers in v at 
NNLO in coupling constant j§|. Thus NRQCD can accommodate these corrections. Hence 
our main aim in this paper is to prove that eq. ((ID is valid at all order in coupling constant. 

Note that in NRQCD an ultraviolet cutoff A ~ M is introduced jl]. Hence the ultraviolet 
(UV) behavior of QCD and NRQCD differ. However, the infrared (IR) behavior of QCD 
and NRQCD remains same 11]. Hence the infrared behavior in NRQCD can be obtained by 
studying the corresponding infrared behavior in QCD. Since the matrix element of the type 
< x}K n £ ( c^ H a hK' n x > is the non-perturbative NRQCD matrix element, it is natural to 
study its infrared behavior at all order in coupling constant by using path integral method. 
Hence we will use path integral method of QCD in this paper. 


IV. HEAVY QUARKS AND THE PATH INTEGRAL FORMULATION OF QCD 

The generating functional in QCD including the heavy quark is given by HQ 


Z[j, r } u , 7] u , rj d , rj d , rj s , r ] s , r ] h , rj h ] = 


] [dip 1 ] [dip 2 ] [dip 2 ] [dip 3 ] [dip 3 ] [c/T] [dV] det( 


d{d^ a ), 

5uj b 


£ i f d 4 x[-±F a l„{Q]-±(d„Q» a ) 2 +J-Q+Y J 3 l=1 [M^d^-mi+gT a 1 »Q“W l +ii l ^ l +ip l r ll \+^i^d„-M+gT a 1 »Q“]*+n h y+'l’Vh] 

( 12 ) 
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where Q^ a is the quantum gluon field, the symbols l — 1, 2, 3 = u, d, s stand for three light 
quarks u, d, s and the symbol h stands for heavy quark and 


KM = W(x) - d„Q‘(x) + gf QAx)Ql(x), 


F a lM = e“„[Q]F'“'“[Q]. 


(13) 

In eq. (1T2|) the fj u , fjd , fj s are external sources for u, d, s quark fields respectively and fjh is 
the external source for the heavy quark field and the term S(d ^\ is the derivative of the 
gauge fixing term under an infinitesimal gauge transformation 


10 , 


12 ] 


= gf abc Q^ b uj c + dL l LU a . 


(14) 


Note that the determinant det( ^%^ - ) in eq. (fl2|) can be expressed in terms of path 
integration over the ghost fields 10]. However, we will directly work with the determinant 
det («MP) in eq. O. 


For the heavy quark Dirac field 4/(cc), the non-perturbative matrix element of the type 
< 0|'F(a:)O n d / (x)'F(a: , )(9( l 'k(a: , )|0 > in QCD is given by 8] 

< 0|'F(x)O.„^(x)^(x / )O;^(x / )|0 > 


— 

5r] h (x) n Srj h (x) 5r] h (x') n 5r] h (x' 


ZfF, ( nr'\ Vui Vui Vdi Vd .1 Vs, Vs, Vh, Vh\ | J=r/u =Vu=Vd =Vd=Vs =Vs =Vh =Vh =0 

(15) 


if the factors O n and 0' n are independent of quantum fields where the suppression of the 
normalization factor Z[ 0] is understood as it will cancel in the final result (see eq. (1133jl ). 


V. INFRARED DIVERGENCES IN NRQCD AND LIGHT-LIKE WILSON LINE 
IN QCD 

The gauge transformation of the quark field in QCD is given by 

^(x) =e i9Ta “ a Wi/;(x). (16) 

Hence one finds that the issue of gauge invariance and factorization of infrared divergences 
in QCD can be simultaneously explained if oj a (x) can be related to the gluon field H Ma (x). 
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Before proceeding to the issue of gauge invariance and the factorization of infrared diver¬ 
gences in QCD let us first discuss the corresponding situation in QED. The gauge transfor¬ 
mation of the Dirac held of the electron in QED is given by 


ip'(x) = e““ (1 >v>(z). 


(17) 


Hence we can expect to address the issue of gauge invariance and factorization of infrared 
divergences in QED simultaneously if we can relate the uj(x) to the photon held A^(x). 

In QED the infrared (or soft) divergence arises only from the emission of a photon for 
which all components of the four-momentum are small. The Eikonal propagator times the 
Eikonal vertex for a soft photon with momentum k interacting with a light-like electron 


moving with four momentum is given by 


fl 9, 


13 


21 ] 


jf 


V 1 


= e 


(18) 


p ■ k + ie l ■ k + ie 

where is the four-velocity of the light-like electron. Note that when we say the ” light-like 
electron” we mean the electron that is traveling at its highest speed which is arbitrarily close 
to the speed of light (|Z| ~ 1) as it can not travel exactly at speed of light (|/j = 1) because 
it has finite mass even if the mass of the electron is very small. From eq. (fTBj) we find 

r d A k l ■ A(k) . f°° „ r d A k 


(27 t) a l ■ k + ie Jo J (27r) 4 


e il ' kX l ■ A{k ) = ie / dXl ■ A(IX) 


/o J Jo 

where the photon held A^(x) and its Fourier transform A^{k) are related by 

d A k 


(19) 


A^(x) = 


(27r) 


-A^(k)e 


ik-x 


( 20 ) 


From eq. (ITDl) we find 


roo r 

ie J dXl ■ A(IX) = i J d A xJ^(x)A ll (x) 


( 21 ) 


where the eikonal current density J A ‘(x) for the light-like charge e is given by 


J*{x) = eV l 


x — IX). 


( 22 ) 


Now consider the corresponding Feynman diagram for the infrared divergences in QED 
due to exchange of two soft-photons of four-momenta k± and k%. The corresponding Eikonal 
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contribution due to two soft-photons exchange is analogously given by 


d 4 h d A k 2 


l ■ A{k 2 )l ■ A(k! 


J (27t) 4 (27r) 4 (/ • (ki + k 2 ) + ie)(l ■ k x + ie) 

POO POO 

e 2 i 2 / d \ 2 / dXil ■ A(lX 2 )l ■ AUXi) 

Jo J \ 2 


e i r c 

Jo 


poo 

dX 2 dXil ■ A(lX 2 )l ■ A(lXi). 
Jo 


( 23 ) 


Extending this calculation up to infinite number of soft-photons we find that the Eikonal 
contribution for the infrared divergences due to soft photons exchange with the light-like 
electron in QED is given by the exponential 


je d\l-A(l\) 


(24) 


where F is the light-like four velocity of the electron. The Wilson line in QED is given by 

(25) 


ie f x f dx^Aa(x) 

£ Jx i 


When A ll {x) = A^(IX) as in eq. (1241) then one finds from eq. (T25l) that the light-like Wilson 
line in QED for infrared divergences is given by 22] 


ie f x dx^Au(x) _ —ie f°°dXl-A(x-\-lX)ie [°° dXl-A(lX) 

3 jo — £ jo £ jo 


(26) 


Note that a light-like electron traveling with light-like four-velocity F produces U(l) pure 
gauge potential A^(x) at all the time-space position x 11 except at the position x perpendicular 


to the direction of motion of the electron (l-x — 0) at the time of closest approach 14, 23, 24], 
When A^ix) = A ll {Xl) as in eq. (|24l) we find l ■ x — XI ■ l — X ^ 0 which implies that the 
light-like Wilson line finds the photon held A^(x) in eq. (l24l) as the U(l) pure gauge. The 
U(l) pure gauge is given by 


A^(x) = d^oj(x) 

which gives from eq. (1261) the light-like Wilson line in QED for infrared divergences 


(27) 


a ieuj(x) — ieuj(O) _ ie J* dx^A^x) _ p —ie d\l-A(x+l\) ie J 0 °° d\l-A(lX) 


(28) 


which depends only on end points 0 and x M but is independent of the path. The path 
independence can also be found from Stokes theorem because for pure gauge 

F^(x) = & l A v {x) - ff'A^x) = 0 (29) 
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which gives from Stokes theorem 


e ief c clxf‘A /l (x) _ e ie f g dy^dx"F^{x) _ ^ 


(30) 


where C is a closed path and S is the surface enclosing C. Now considering two different 
paths L and M with common end points 0 and x M we find 


e ie § c dxtlA A x ) — e ie S L dxlJ 'A )1 (x)—ie J M dx^A^x) _ ^ 


(31) 


which implies that 


e ie Jo dxl * A A x ) 


(32) 


depends only on end points 0 and x M but is independent of path which can also be seen from 
eq. (j28jl . Hence from eq. (1281) we hnd that the abelian phase or the gauge link in QED is 
given by 


ie dAl-A(x+lA) _ gieoj(x) 


(33) 


From eqs. (fI71) and (j33l) one expects that the gauge invariance and factorization of infrared 
divergences in QED can be explained simultaneously. 

One can recall that the gauge invariant greens function in QED 


r x 2 


G(xi, X2) —< 0\ij}(x2) x exp[ie / dx^A^x)] x ijj(xi)\0 > 


(34) 


' X\ 


in the presence of background field A^(x) was formulated by Schwinger long time ago 25]. 


When this background held A IJ, (x) is replaced by the U(l) pure gauge background held as 
given by eq. (127|) then one finds by using the path integral method of QED that 


e ie^ X2 ) < 01 ^) ^( Xl )|0 > A e~ ieu(xi) =< 0$(x 2 ) ?p(xi)\0 > 

= e -ief 0 °°dAl-A(x 2+ lA) < ^ Xl )\0 >A e ief 0 °°dAl A {xi+ lA) ( 35 ) 


which proves the gauge invariance and factorization of infrared divergences in QED simul¬ 
taneously. I 11 eq. (I35|) the < O^a^) > is the full Green’s function in QED and 

< O|'0(m 2 ) ^(aq)|0 >a is the corresponding Green’s function in the background field method 


of QED. This path integral technique is also used in 
divergences in non-equilibrium QED. 


18] to prove factorization of infrared 
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Hence we find that the gauge invariance and factorization of infrared divergences in QED 
can be studied by using path integral method of QED in the presence of U(l) pure gauge 
background held. Therefore one expects that the gauge invariance and factorization of 
infrared divergences in QCD can be studied by using path integral method of QCD in the 
presence of SU(3) pure gauge background held. 

Now let us proceed to QCD. In QCD the infrared (or soft) divergence arises only from 
the emission of a gluon for which all components of the four-momentum are small. The 
Eikonal propagator times the Eikonal vertex for a soft gluon with momentum k interacting 


with a light-like quark moving with four momentum p M is given by 


13-19 




gT° 


p*- 


— gT a 


V 1 


(36) 


p ■ k + ie l ■ k + ie 

where P is the four-velocity of the light-like quark. Note that when we say the ” light-like 
quark” we mean the quark that is traveling at its highest speed which is arbitrarily close to 
the speed of light (|Z| ~ 1) as it can not travel exactly at speed of light (|Z| = 1) because 
it has finite mass even if the mass of the light quark is very small. On the other hand the 
gluon is massless and hence it always travels at speed of light and is exactly light-like. From 
eq. (136|) we find 


gT a 


d A k l ■ A a (k) 
(27r) 4 l ■ k + ie 


= -gT a i 


dX 


d 4 k 

Wr 


e“"I • 4” (In = igT a / d\l ■ A‘(l A) 


where the gluon held A^ a (x) and its Fourier transform A^ a (k) are related by 

d A k 


A^ix) = 


(2tt) 


; A^(k )e 


ik-x 


(37) 


(38) 


Note that a path ordering in QCD is required which can be seen as follows, see also 20]. 


The Eikonal contribution for the infrared divergence in QCD arising from a single soft- 
gluon exchange in Feynman diagram is given by eq. (I37|h Now consider the corresponding 
Feynman diagram for the infrared divergences in QCD due to exchange of two soft-gluons 
of four-momenta k\ and k%. The corresponding Eikonal contribution due to two soft-gluons 
exchange is analogously given by 

2 /• d A k i d 4 k 2 T a l ■ A a (k 2 )T b l ■ H fe (Aq) 


g 


J (27r) 4 (27 t) 4 (l ■ (A* + k 2 ) + ie)(l ■ + ie) 

n oo poo 

g 2 i 2 / d\ 2 / dX\T a l ■ A a (lX 2 )T b l ■ A b (l AQ 

Jo J a 2 

-9 2 p oo poo 

—V / dX 2 / dX\T a l ■ A a (lX 2 )T h l ■ A b (l AD 

! Jo Jo 


9 

2 ! 


( 39 ) 
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where V is the path ordering. Extending this calculation up to infinite number of soft- 
gluons we find that the Eikonal contribution for the infrared divergences due to soft gluons 
exchange with the light-like quark in QCD is given by the path ordered exponential 


/»oo 

V exp [ig J dXl ■ A a {lX)T a 


(40) 


where is the light-like four velocity of the quark. The Wilson line in QCD is given by 


Ve 


ig /7 dx»A“(x)T a 


which is the solution of the equation 26] 


d^Six) = igT a A a (x)S(x) 


(41) 


(42) 


with initial condition 


S(xi) = 1 . 


(43) 


When A txa (x) = A^ a (l\) as in eq. (jjOj) we find from eq. (ITT]) that the light-like Wilson line 


in QCD for infrared divergences is given by 


22 ] 


TV 9 Jo *dx»A°(x)T“ = \p e -i9 f 0 °° dM-A“(x+l\)T a 1 jjjg f 0 °° d\l-A b (l\)T 


(44) 


A light-like quark traveling with light-like four-velocity produces SU(3) pure gauge 

potential A IM (x) at all the time-space position x M except at the position x perpendicular to 


14 


24] 


the direction of motion of the quark (l ■ x = 0) at the time of closest approach 
When A /J,a (x) = A ga (Xl) as in eq. (ITU]) we find l ■ x = XI ■ l = X ^ 0 which implies that the 


light-like Wilson line finds the gluon field A^ a (x) in eq. (ITU]) as the SU(3) pure gauge. The 
SU(3) pure gauge is given by 


T a A“(x) = -[d,U(x)} U~\x), 
ig 


U(x) = e igTaujaix) 


which gives 


U(x f ) = Ve a G d ’“ Ai( ’ )T "u(x t ) = 


(45) 


(46) 


Hence when A /J,a (x) = A g,a (Xl ) as in eq. (]TU1) we find from eqs. (ITTj) and (ITUl) that the 
light-like Wilson line in QCD for infrared divergences is given by 

Ve ig^dx^A-{x)T- = e igT-u-{x) e -igT b oj\ 0 ) = p e ~ig /“ dM-A“(x+l\)T° pjg f Q °° d\l-A b {l\)T b 
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which depends only on end points 0 and x M but is independent of the path. The path 
independence can also be found from the non-abelian Stokes theorem which can be seen as 
follows. The SU(3) pure gauge in eq. (1451) gives 


K„\A] = d„Al(x) - 8Mx) + gf^A b Jx)AUx) = 0. 


(48) 


Note that from eq. (1481) we hnd the vanishing physical gauge invariant field strength square 
F liua [^ 4 ] [H] w hen A tm (x) is the SU(3) pure gauge as given by eq. (1451) . Hence in classical 
mechanics the SU(3) pure gauge potential does not have an effect on color charged particle 
and one expects the effect of exchange of soft gluons to simply vanish. However, in quantum 
mechanics the situation is a little more complicated, because the gauge potential does have an 
effect on color charged particle even if it is SU(3) pure gauge potential and hence one should 
not expect the effect of exchange of soft gluons to simply vanish [bf]. This can be verified by 
studying the non-perturbative matrix element in QCD such as < T(a;)T(a: / )'^ , (x , ') v I / (x ,,, )... > 
in the presence of SU(3) pure gauge background field. 


Using eq. (148|) in the non-abelian Stokes theorem 27] we hnd 

p e W J* dx ' X A\(x')T b 


Ve i 9 § c dx»A«{ x )T° =Vexp [ ig j dx P dx - 


KMT a 


Ve i9 -C dx " s A c s (x")T c 


= 1 


(49) 


where C is a closed path and S is the surface enclosing C. Now considering two different 
paths L and M with common end points 0 and x M we hnd from eq. (j49li 


Ve ig§ c <WA%{x)T* = 'Pexplig J^ dx »A“(x)T a - ig j dx^A a ^(x)T a ] 

= [Ve i9 $L 


dx 11 A^(x)T a 


V e - i9 i M dxVAb AAT b 


(50) 


which implies that the light-like Wilson line in QCD 

Ve ig f* dx^A-(x)T a 


(51) 


depends only on the end points 0 and x M but is independent of the path which can also be 
seen from eq. (1471) . Hence from eq. (FT7D we hnd that the non-abelian phase or the gauge 
link in QCD is given by 


= p e -i9 /“ d\l-A“(x+l\)T- = e igT°u“(x)' 


(52) 
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In the adjoint representation of SU(3) the corresponding path ordered exponential is given 


by 


Vexp{-ig r dXl ■ A c (x + IX)T (A)c } = e ^ T(A)c " c (*) ; ( T (A)c ) ab = -if abc . 

Jo 


(53) 


To summarize this, we find that the infrared divergences in the perturbative Feynman 
diagrams due to soft-gluons interaction with the light-like Wilson line in QCD is given by 
the path ordered exponential in eq. (1401) which is nothing but the non-abclian phase or the 
gauge link in QCD as given by eq. (I52|) where the gluon held A^ a (x) is the SU(3) pure gauge, 
see eqs. (l45]h (l46]h fl47|h This implies that the effect of soft-gluons interaction between the 
partons and the light-like Wilson line in QCD can be studied by putting the partons in the 
SU(3) pure gauge background held. Hence we find that the infrared behavior of the non- 
perturbative matrix element such as < 0|T(x)T(x / )'l , (x")'F(x ,, ')...|0 > in QCD due to the 
presence of light-like Wilson line in QCD can be studied by using the path integral method 
of the QCD in the presence of SU(3) pure gauge background held. 


It can be mentioned here that in soft collinear effective theory (SCET) 28] it is also 

“ y t0 use the r M f a of background flelds H wdl drftaed — to SCTeral 

distinct gluon helds 1151 1. 

n 

As mentioned earlier, in NRQCD an ultraviolet cutoff A ~ M is introduced pj. Hence the 
ultraviolet (UV) behavior of QCD and NRQCD differ. However, the infrared (IR) behavior of 


QCD and NRQCD remains same 


11| . Hence the infrared behavior in NRQCD can be stud¬ 


ied by studying the corresponding infrared behavior in QCD. Hence we find that the infrared 
behavior of the non-perturbative NRQCD matrix element < Q\x^K n £(a' H a,H)^K' n x\® > hi 
eq. ([I]) can be obtained by studying the infrared behavior of the non-perturbative matrix 
element in QCD of the type < 0|T(a;)O n T(a;)T(a; , )O, , l v h(a; , )|0 > where O n , 0' n are appropri¬ 
ate factors which identify the state of the heavy quark-antiquark system such as the color 
singlet state or color octet state etc.. 

Note that a massive color source traveling at speed much less than speed of light can 
not produce SU(3) pure gauge held 14j, [23|, [24]. Hence when one replaces light-like Wilson 
line with massive Wilson line one expects the factorization of infrared divergences to break 
down. This is in confirmation with the finding in 29] which used the diagrammatic method 
of QCD. In case of massive Wilson line in QCD the color transfer occurs and the factorization 
breaks down. 


14 











VI. EIKONAL CURRENT OF THE LIGHT-LIKE CHARGE GENERATES PURE 
GAUGE FIELD IN QUANTUM FIELD THEORY 


In order to study factorization of infrared divergences by using the background field 
method of QED, the soft photon cloud traversed by the electron is represented by the pure 
gauge background field A^(x) js] due to the presence of the light-like Wilson line, where one 
represents the quantum photon field by Q^(x). As mentioned above, in classical mechanics 
the assertion that the gauge field that is produced by a highly relativistic (light-like) particle 


is a pure gauge 


14, [23j, 24], One may ask a question if this assertion is correct in quantum 


field theory. In this section we will show that this assertion is correct in quantum held 
theory. We will use path integral formulation of the quantum held theory for this purpose. 

The generating functional for the gauge held in the quantum held theory in the presence 
of external source J M (x) in the path integral formulation is given by 


m = 

where is the quantum photon held and 

F^[Q]=d^Q u (x)-d u Q^x), 


The effective action S e ff[J] is given by [34] 


< 0|0 >j= 


ZjJ] 

Z[ 0 ] 


KM = F^[Q]F^[Q], 


— g ®^e// U] 


where 


S e ff[J\ = — -J d A xd A x' J^(x)D^ lv (x — x') J v (x') 


where D lw (x — x') is the photon propagator. 

The photon propagator in the coordinate space is given by 

D^(x - x') = -^[g^u + ^ d^d u \5 (/L) (x - x'). 


Using eq. (T58]) in (157]) we find 

Seff[J] = ~\f d 4 xJ^(x)-^[g^ + 

From the continuity equation we have 


(a ~ !) 
& 2 




(55) 


(56) 


(57) 


(58) 


(59) 


d r P(x) = 0. 


(60) 
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(61) 


Using eq. (1601) in (|59|) we find 


Seff[J] = ~\J d A xJ^(x)^J^x) 


A. Derivation of Coulomb’s Law For Static Charge in Quantum Field Theory 

First of all, by using the path integral formulation of the quantum field theory we will 
derive Coulomb’s law for static charge. Note that the derivation of the Coulomb’s law 
by using path integral formulation of the quantum field theory is not necessary to prove 


factorization theorem. We have included it here only to demonstrate the correctness of the 


prediction of the path integral formulation in quantum field theory which we will use (see 
below) to show that the eikonal current of the light-like charge generates pure gauge field in 
quantum held theory. 

In order to derive Coulomb’s law by using path integral formulation of the quantum 


held theory we consider two static charges at positions X and X' respectively. The current 
density for this two static charges is given by 


J“{x) = eJ“°<5 <3 >(x -X) + eS^S<- 3 >(x - X'). 


(62) 


Using eq. ()62l) in (T6T1) and neglecting the self energies we find in the time interval t that 


S?MA = e 1 [dtX 6 ( 3)(X _ X') = -t e ; AvU 1 ) = (63) 


which gives the (effective) potential energy V e ff[J] of the interaction between two static 
charges to be 



(64) 


which reproduces the Coulomb’s law. Hence we have shown that the assertion that a charge 
at rest generates a Coulomb gauge held is correct in quantum held theory. 

B. Effective Lagrangian Density of Light-Like Eikonal Current in Quantum Field 
Theory 

Similarly using the above procedure in quantum held theory we will show that the as¬ 
sertion that a light-like charge generates pure gauge held is correct in quantum held theory. 
This can be shown as follows. 
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The eikonal current density of the charge e with light-like four-velocity P is given by eq. 
By using the path integral formulation of the quantum field theory we find by using 

eq. HT2b in (IMP that for light-like eikonal current the effective lagrangian density is given by 

.2 (/2) 2 




2 (l ■ x) 4 ' 


For light-like four-velocity we have 


l 2 = PL = 0 . 


(65) 


( 66 ) 


Hence from (j65j) and 
density is given by 


we find that for light-like eikonal current the effective lagrangian 


£e//0) = 


l ■ x 7 ^ 0 


(67) 


at all the time-space position x M except at the spatial position perpendicular to the motion 
of the charge (l ■ x — 0 ) at the time of closest approach (x 0 = 0 ). 


C. Interaction Between Non-Eikonal Current and the Gauge Field Generated by 
Light-Like Eikonal Current in Quantum Field Theory 

Similarly by using the above path integral formulation calculation we find from eq. 


that the interaction between the (light-like or non-light-like) non-eikonal current and the 
gauge field generated by the light-like eikonal current gives the effective (interaction) la¬ 
grangian density 


4//(*0 = i 2 


e 2 (l ■ q)(q ■ x) — {l ■ x)q 2 


( 68 ) 

2 (l ■ x) 3 [(q ■ x) 2 — q 2 x 2 } 2 
where is the (light-like or non-light-like) four-momentum of non-eikonal current of charge 
e and P is the light-like four-velocity of the eikonal current of charge e. 

For light-like eikonal current we find from eqs. (1661) and (1681) that effective (interaction) 
lagrangian density due to the interaction between the (light-like or non-light-like) non-eikonal 
current of four-momentum and the gauge field generated by the light-like eikonal current 
of four-velocity P is given by 


C:/ f (x) = 0 , for 

This is also obvious from eq. (f 8 T|i . 


q ■ x 7 ^ 0 , 


l ■ x 7 ^ 0 . 


(69) 
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D. Pure Gauge Field Generated By Eikonal Current of Light-Like Charge in 


Quantum Field Theory 


Hence from eqs. (1571) and (155|) we find that the eikonal current for light-like charge 
generates pure gauge field in quantum field theory. From eqs. (157|) and (155|) we find that 


the assertion that a light-like charge generates a pure gauge field is correct in quantum 
theory which is consistent with the corresponding result in classical mechanics 14, 23], 


he 


24], 


VII. PURE GAUGE FIELD IN QUANTUM FIELD THEORY DESCRIBES SOFT 
(INFRARED) DIVERGENCE 


In this section we will show how the pure gauge field is used in quantum field theory to 
describe soft (infrared) divergences. Consider an incoming electron of four momentum q M 
and mass m emitting a real photon of four momentum k C The corresponding Feynman 


diagram contribution is given by 

1 


30] 


M = 


where we write 


r ) v q v — — m 


7 ^(k)u( g ) = - q -^u {q) + k2^m u(9) (70 ) 


q ■ k 


2 q ■ k 


and 


. . q-e(k) 

eikonal =- JTu{q) 

Q * /C 


-Adnon—eikonal — ^ , 'U J \Q) • 


(71) 


2 q ■ k 


(72) 


From eq. (3.2) of [3Q1] we write the gauge field as 


<?(k) = khh - = ej) hys (fc) + t» n Jk) 


where 


=(fe) = KW - v 


Q-e(.k) 

q ■ k 


(73) 


(74) 


is the physical gauge field [corresponding to transverse polarization of the gauge field] and 

,q-e{k) 


^ P ur e (k) = k»- 


q ■ k 


( 75 ) 


is the pure gauge field [corresponding to longitudinal polarization of the gauge field]. 
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Now using eq. (J73j) in eq. (1701) we find that the total contribution of the Feynman diagram 
is given by 


M 


M 


eikonal 


+ M 


non—eikonal 


(76) 


where 


« y Q ' e phys(fc) / \ Q ' e pure(^) / \ Q ' e pure(^) / n. 

Alcikonal =-7-«(<?)--«(?) =-^- U{q) 


and 


M non—ei 


non—eikonal 


q ■ k 

fc l/ 7^^h ys ( fc ) 


q ■ k 


q ■ k 


(77) 


, * , V7„7 ^ (k) fc"7„7(fc) 

2, • k U(q) + -Vi- U{q) =-2?^- U(q) ■ (78) 

Hence in the soft photon limit (k 0 , ki, k 2 , k 3 ) 4 0 we find from the eqs. (1701) and (1771) 
that 


-M 


eikonal 


<?-e P ur e(fc) 


q ■ k 


-u(q) = 


q ■ k 


u(q) -4- oo 


as (fco> fcij ^ 2 ) k%) —> 0(79) 


which implies that the physical gauge held [corresponding to transverse polarization] does 
not contribute to the soft (infrared) divergences in quantum held theory and the soft (in¬ 
frared) divergences can be calculated by using pure gauge held [corresponding to longitudinal 
polarization] in quantum held theory. 

Similarly from eqs. (1701) and (17H1) we hnd that 


-^drion—' 


non—eikonal 


= kU ^yp^ k) u(q) -4 finite 


2 q ■ k 


2 q ■ k 


as 


(k 0 , k\, &2, k 3 ) —> 0 


(80) 


which contribute to the finite part of the cross section which implies that pure gauge held 
[corresponding to longitudinal polarization] does not contribute to the hnite cross section 
and the hnite cross section can be calculated by using physical gauge held [corresponding 
to transverse polarization]. 

Hence we hnd that the non-eikonal-line part of the diagram as given by eq. (1801) is 
necessary if we are calculating the hnite value of the cross section but it is not necessary if 
we are calculating the relevant infrared divergence behavior. The relevant infrared divergence 
behavior can be calculated by using the eikonal-linc part of the diagram as given by eq. (1791) . 

For this reason, in the proof of NRQCD factorization of infrared divergences for heavy 
quarkonium production at NNLO in coupling constant, the non-eikonal-line part of the 
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diagram as given by eq. (ISOl) is not considered as the full calculation of the cross section or 
fragmentation function at NNLO will be daunting but fortunately the analysis of relevant 
infrared behavior at NNLO requires only the eikonal approximation as given by eq. (1791) . 
see the discussion in the last paragraph of section 4 of 4]. Similarly the full calculation of 
the cross section or fragmentation function at all order in coupling constant by using the 
non-eikonal-line part of the diagram as given by eq. (180]) will require non-perturbative QCD 
which is not solved yet but fortunately the analysis of relevant infrared behavior at at all 
order in coupling constant requires only the eikonal approximation as given by eq. d79|) . 

Hence we find that we do not need to calculate the finite value of the cross section (or the 
full cross section) [which will require the non-eikonal-line part of the diagram as given by eq. 
(189]) ] to study the relevant infrared divergence behavior. The relevant infrared divergence 
behavior can be calculated by using eikonal approximation as given by eq. (179]) . 

From eq. (180]) we find that 


M 


pure gauge field _ ^ 7i'7/Rpure(^0 


non—eikonal 


2 q ■ k 


u(q ) = 0. 


(81) 


We are interested in the infrared divergence behavior due to the presence of the light-like 
Wilson line. We have shown in eqs. (167]) and (169]) that the eikonal current of the light-like 
charge generates pure gauge held in quantum held theory. Hence from eqs. (l67j) . (1791) . 
(]69]) and (TSTT) we find that the soft (infrared) divergence behavior due to the presence light¬ 
like Wilson line can be studied by using pure gauge held in quantum held theory without 
modifying the finite value of the cross section. 


VIII. HEAVY QUARK-ANTIQUARK NON-PERTURBATIVE MATRIX ELE¬ 
MENT IN THE PRESENCE OF LIGHT-LIKE WILSON LINE IN QCD 

We have seen in section V that the infrared behavior of the non-perturbative NRQCD 
matrix element < 0|yf K n fi(a} H aH)t}K' n x\0 > in eq. ([I]) can be obtained by study¬ 
ing the infrared behavior of the non-perturbative matrix element in QCD of the type 
< 0|T(a;)(9 n T(x) 1 F(a; / )O( l T(a; , )|0 > where O n , 0' n are appropriate factors which identify 
the state of the heavy quark-antiquark system such as the color singlet state or color octet 
state etc.. Similarly, we have also seen in section V that the infrared behavior of the non- 
perturbative matrix element in QCD of the type < 0|4 , (x)0„'I / (a;) 1 F(x / )0( l T(x / )|0 > due 
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to the presence of light-like Wilson line in QCD can be studied by using the path integral 
method of the QCD in the presence of SU(3) pure gauge background field. Hence we use the 
path integral formulation of the background field method of QCD to study non-perturbative 
matrix element < 0|'k(x)0„ 1 k(a;)'l , (x')0^ v h(x')|0 > in QCD in the presence of SU(3) pure 
gauge background field as given by eq. (PX3T) . 


Background field method of QCD was originally formulated by ’t Hooft 31] and later 


extended by Klueberg-Stern and Zuber 


32, 


33( ] and by Abbott k2j. This is an elegant for¬ 


malism which can be useful to construct gauge invariant non-perturbative green’s functions 


in QCD. This formalism is also useful to study 


quark and gluon production from c 


assical 


chromo field 34] via Schwinger mechanism 35h to compute /3 function in QCD 36], to 
perform calculations in lattice gauge theories |37] and to study evolution of QCD coupling 


constant in the presence of chromofield 


38]. 


In the background field method of QCD the generating functional is given by Q,H, 

,6G a (Q) 


32] 


Z[A, J, t), rj\ — J [dQ][di/j][di/j] det(- 


5u b 


-) 


a i f d 4 xl-^F^lA+Q]-^(G a (Q)) 2 +^a li -m+ 9 T a 7 f‘(A+Q)^+J-Q+f!^r ) ] 


(82) 


where Q^ a (x) is the quantum gluon field and the gauge fixing term is given by 

G a (Q) = d.Cy* + gf abc A\Q» c = D /1 [A}Q lia (83) 

which depends on the background field A ,,a and 

KM + Q\ = 8M: + «] - 8AK + «] + sKM* + Ql ]K + «]• (84) 


We have followed the notations of 


0,0,E| 


and accordingly we have denoted the quantum 


gluon field by Q^ a and the background held by A^ a . The determinant c 


et(^) in eq. 


32]. However, we will 


can be expressed in terms of path integration over the ghost fields 10, 
directly work with the determinant det( ; ^^) in eq. (JS2D- 

Note that the gauge fixing term A^(G a (Q))' 2 in eq. (j82l) [where G a (Q ) is given by eq. 
is invariant for gauge transformation of A“: 


<5A“ = gf abc A b ^u> c + d ll oj a , (type I transformation) 


(85) 


provided one also performs a homogeneous transformation of Q 


; Q& 


s<z = ?/*&•. 


( 86 ) 


21 






































The gauge transformation of background field A a as given by eq. ([85]) along with the 
homogeneous transformation of in eq. (t 86 |) gives 

KK + Q“) = gr bc (A^ + Q \> c + d,u a (87) 

which leaves • m \F a2 flv [A + Q\ invariant in eq. (j82j) . 

For fixed i.e., for 


8A a = 0 , (type II transformation) 


the gauge transformation of 



SQl = gf abc (Al + Q> c + d^us 


( 88 ) 


(89) 


gives eq. (1ST]) which leaves —\F a2 [A + Q] invariant in eq. (]82]h 

Extending eq. (|82|) to include heavy quark [by using the lagrangian density from eq. (J6])] 
we find that the generating functional in the background field method of QCD is given by 

f - - - 8G a (Q ) 

Z[A, J , rj u , fj u , rj d , fj d , r} s , fj s , rj h , fj h \ = J [dQ} [dipi] [di>i\ [dip 2 ] [# 2 ] [dip 3 ] [dip 3 \ [dT] [d'F] det( ^ • ) 

exp[ S j + Q}~ ^(G“(Q )) 2 + J-Q 

+ - m + gT a ^(A + QY^ipi + fhipi + 

i=l 

+T[* 7 ^ - M + + Q)“]T + 77 ^ + T % ]]. (90) 


Note that in the absence of external sources a pure gauge can be gauged away from the 
generating functional. However, in the presence of external sources a pure gauge can not be 
gauged away from the generating functional. It is useful to remember that, unlike QED jsi], 
finding an exact relation between the generating functional Z[J,rj u ,fj u ,rjd,fjd,r} s ,fj s ,rjh,rih\ 
in QCD in eq. (TT21) and the generating functional Z[A, J, rj u , fj u , rj d , fjd, Vs, Vs, Vh, Vh] in the 
background field method of QCD in eq. (jUUl) in the presence of SU(3) pure gauge background 
field is not easy. The main difficulty is due to the gauge fixing terms which are different in 
both the cases. While the Lorentz (covariant) gauge fixing term —■^-(d lx Q^ a ) 2 in eq. (TT 2 l) 
in QCD is independent of the background field A^ a (x), the background field gauge fixing 
term — ^(G a (Q)) 2 in eq. (1901) in the background held method of QCD depends on the 


background held H M “(x) where G a (Q) is given by eq. 


[i3, El, 


32 


Hence in order 
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to study non-perturbative matrix clement in the background field method of QCD in the 
presence of SU(3) pure gauge background field we proceed as follows. 

By changing Q —> Q — A in eq. (15U1) we find that 


Z[A, J, rj u , r) u , rj d , g dl rj s , r} a , rj h , rj h ] 

= e ~ l f d4xJ ' A J [dQ][d^i][d^i][d^][# 2 ][#3][#3][^][d^] det ( ~^~^ ) 

e < / d A x[-\F-l v [Q]-^{G a s {Q)) 2 +J-Q+Y^ =1 [Mi^d^- mi +gT a ^Ql}i )l +n l ^ l +m^+nh^d^-M+gT^Q^+r lh <i!+^i lh } 

(91) 


where the gauge fixing term from eq. (t83|) becomes 


G){Q) = d lt Q IJa + gf abc A b ^ c - d,A ,M = D^AlCy 0 ' - d^, (92) 

and eq. (j86|) [by using eq. (185]) . type 1 transformation jla Idl becomes 


SQl = gf^QW + d„<S. 


(93) 


The eqs. (192)1 and (1991) can also be derived by using type II transformation which can be 
seen as follows. By changing Q —> Q — A in eq. (190)1 we find eq. (191)1 where the gauge fixing 
term from eq. (183)1 becomes eq. Q92j) and eq. (189|) [by using eq. f[88j) ] becomes eq. f)93]) . 
Hence we obtain eqs. (ED), ([92]) and (l93|) whether we use the type I transformation or type 
II transformation. Hence we find that we will obtain the same eq. (1120)1 whether we use the 
type I transformation or type II transformation. 

Note that 


A!‘(x) = A‘{x) + gf^V(x)A%x) + B^‘(x) 


(94) 


in eq. (185)1 is valid for infinitesimal transformation (c j « 1) which is obtained from the 
finite equation 

T a A'“(x) = U(x)T a A“(x)U-\x) + —[dyj(x)}U~ l (x) 1 U(x) = e < ® T ° w ° (a!) . (95) 

ig 

Simplifying infinite numbers of non-commuting terms we find 


e ~igT b u b (x) rpa ^,gT c u c {x) 


ij 



where 


(96) 


M ab (x) = f^u^x). 


(97) 
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23J we find that 


Hence from eqs. (l95]h (196|) and 

p gM{x) _ i 

KM = + I gM(x) U IW) 1 (98) 

where M ab (x) is given by eq. (I97|) . Similarly, the equation 

O'; ?(*) = + 9 f abc oo\x)Q b ^x) + VC*0 (99) 

in eq. (193|) is valid for infinitesimal transformation (a; << 1) which is obtained from the 
finite equation 

T a Q’“(x ) = U(x)T a Ql(x)U-\x) + -{d^U{x))U~\x) (100) 

ig 

which gives 

p gM(x) _ -| 

<r»w = [e 9 " (i) u<?fw + i gM(x) u (Vmi (ioi) 

where M ab (x) is given by eq. dHZD- 

Changing the variables of integration from unprimed to primed variables in eq. (19TT) we 
find 


Z[A, J, rj u , g Ul g d , rj d , rj s , t] s , r) h , Vh] 


_ —i J d 4 xJ-A 


'} [dip[\ [dip[] [dip 2 ] [dip' 2 \ [dip' 3 \ [dip' 3 \ [d'h / ] [dT'] det ( 


5G%Q'), 

~s^~- 




( 102 ) 


This is because a change of variables from unprimed to primed variables does not change 
the value of the integration. 

Under the finite transformation, using eq. (11011) . we find 

r)(l ,a 

[dQ’} = [dQ\ det[_^-] = [dQ\ det[[e 9M(T) ]] = [dQ]exp[Tr(hi[e 9M(x) ])] = [dQ\ (103) 
where we have used (for any matrix H ) 

det# = exp[Tr(lnif)]. (104) 

Similarly the fermion fields transform accordingly, see eq. (USD, he., 

iP[(x) = e^^ip^x), ^'(x) = e i9Ta “ a( - x) V(x). (105) 
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Using eqs. (j 101 [) and (I105j> we find 


= [# 'i][<#i], [# 2 ] [# 2 ] = [# 2 ] [# 2 ], [# 3 ] [# 3 ] = [dV> 3 ][# 3 ], 

[d*W] = [d*][d*], tiiiYd, -rm + gT a ^Q'M = - rrg + gT a YQ a ^u 

*'[* rd, -M + gT a YQ'ZW = - m t + gT a ^Q \*]*, F a2 ^[Q'] = 

(106) 


Using eqs. (I103P and (I106j> in eq. (I102|> we find 


Z[A, J, T) u , T) u , gd, rjd, Vs, Vs, Vh, Vh\ 

= g — * f d 4 xJ-A J [dg] [d^] [d^] [d^ 2 ] [d^; 2 ] [d^ig ][di/r 3 ][d^][d^] det{ SG ^^ ) 

g i J d A x[-\F-l v [Q)-^(G a f {Q l )) 2 +J-Q'+Y^ = 1 [Mi^d^-rn l +gT^Ql]i )l +n l ^[+^[m}+nira^-M+gT^Q^+f lh ^+^ Vh }_ 

(107) 


From eq. (192|) we find 


GJ(Q') = a,/:/"" + gf^A^Q'i* - a„4'"'. 


(108) 


By simplifying the infinite number of non-commuting terms in the SU(3) pure gauge in eq. 
we find 123] 


A» a {x) = d^oj b (x) 


e gM(x) _ ^ 

gM (x) 


(109) 


ab 


where M n &(:r) is given by eq. (J97} . By using eqs. (1101 1) and (11091) in eq. (11081) we find 

p gM(x) _ 1 

G%Q') = d^He^%„Q b Jx) + l , U PXWI] 


+gf abc [d fl uj e (x) 


e gM(x) _ 

gM (x) 


gM(x) 

p gM(x) _ 1 

nie^uotix) + i e u [d^ d m 


be 


-d,[d^ b (x) 


e gM{x) _ ^ 


gM (x) 


( 110 ) 


- ab 


which gives 


G%Q') = SV^UQ^x)] 

'g#W _ 2 


+gf abc [d^ e (x) 


gM (x) 


p gM(x) _ 1 

Ee^UQtix) + [ gM(x) U (111) 
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From eq. (lllip we find 

Gy(Q') = d^e^UQ^x)] + 

which gives 

Gf{Q') = [e^Ud^x) 


e gM(x) _ ^ 

gM (x) 


iir 1 -^*)] (H2) 


6 e 


+Qj(^)^[[e 9M(a:) ] a6 ] + [<W(*) 


From 


23 ] we find 

d ^ e igT a «“(*)]„ = ^[^(. t )] 


' e gM{x) _ ^ 

gM (x) 

e gM{x) _ 2" 

gM{x) 


W^U^UQtWl (H3) 


be 


T? k {e ,3TVt -%j 


ab 


which in the adjoint representation of SU(3) gives (by using T 6 “ = — if abc ) 

' e gM(x) _ 1" 




gM(x) 


g fbac [g-M(i)] 


ccf 


be 


where M ab (x ) is given by eq. (EH). Using eq. (Ill 5D in (11131) we find 

G a f {Q') = [e^Ud^x) 

which gives 

(Gy(Q ')) 2 = (d,Q^(x)) 2 . 

Since for n x n matrices A and B we have 


det(AB) = (detA)(detS) 


we find from eq. (jl 16j> that 


d e t [^«l] = det| d&^, = Aet[m ^ Ml; 


(114) 


(115) 


(116) 


(117) 


(118) 


= det[ 


det[[e^% c ] 
8{d tl Q m {x)) 


det[ 


8{d^Ql(x)) 

8u b 1 


= exp[Tr(ln[e 9M(3:) ])] det[ 


5(d^ a (x))- 

8oj b 


5(jj b 


(119) 


Using eqs. ()1 17|> and (1 119 [) in eq. ()107j> we find 
Z[A, J, rj u , fj u , rj d , fj d , rj s , rj s , rj h , fj h \ 


_ —i J d 4 xJ-A 


J [dipi] [# 2 ] [# 2 ] [# 3 ] [# 3 ] [d^] [d'h] det [ 


8(d fl Q fia (x))- 

6uj b 


J f d 4 X [-jF^lQ]-^(d M Q^)MJ-Q'+^i = il^ili7^d /1 -m l +gT^Q^ l +fj l i,;+^r, l ]+^[i 7 ^-M+gT^Q^+fj h 9'+^' Vh ] 


( 120 ) 
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From eqs. (j!09() and (1 101 1) we find 


Q'» - A‘Jx) = [e» M «UQ» 


( 121 ) 


where M ab (x ) is given by eq. dHZD- 

Note that eqs. (11201) . (I12ip and (fl6|) are valid whether we use type I transformation [see 
eqs. (155]) and (155]) ] or type II transformation [see eqs. (155]) and (155]) ]. 

However, since eq. (155]) is used to study the gauge transformation of the Wilson line in 
QCD, we will use type I transformation [see eqs. (155]) and (156]) ] in the rest of the paper 
which for the finite transformation gives [12, [32] 


f‘(x) = [e gMi - z) Wl{x) 

where M ab (x ) is given by eq. (157]) . From eqs. (]120|) . (1 121 () and (11221) we find 


( 122 ) 


Z[A, J', r] u , rj u , r] d , rj dl r] s , rj s , r] h , r] h \ 


] [dVh] [dip-. 2 ] [dip 2 ] [dips] [dips] [d^] [d^] det [ 


5(d ll Q txa {x)) 
5u b J 


J f d 4 ‘x[-±F a l„[Q]-±(d l _ L Qi J - a ) 2 +J.Q+J2^ =1 [4’l +^[i7^(9 M -M+sT a 7' i Q£]'I'+r7f l 'I''+'I''r/ /l ] 


(123) 


Under the non-abelian gauge transformation the fermion sources transform as [l2, 3^] 


rf h (x) = e l9Tauja{x) ri h (x). 


V' h *' = fj h 'V, *'r/ h = *Tj h 


(124) 


(125) 


Vi(x) = e wTa “ a{x) rii{x), 

From eqs. (11051) and (11241) we find 

rfitpi = m^Pi, i’lV'i = $iVi, 
which gives from eq. (II 23j) 

Z[A, J', rju, fj'u, rid, v'd, Vs, v' s , Vh, Vh] 

J [#J [# 2 ] [# 2 ] [# 3 ] [# 3 ] [<»] [»] del[KMUhl| 

J J d 4 x[-^F a l v [Q]-^(d^ a ) 2 +J-Q+J2 3 l=1 [Mn ll O l ,-rn l +gT^Q'l l ]^ l +ij l i n +^ l r ]l ]+^[i^d l ,-M+gT a 'yi b Q<l l ]^+g h ^+^ h } 


( 126 ) 


Hence from eqs. (11261) and (1151) we find 


Z[J, rju, riu, rid, rid, Vs, Vs, Vh, Vh] = Z[A, J', r/ u , V u , r/ d , r/ d , 77 ', v' s , v'h, Vh] 


(127) 
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when the background field 2D a (x) is the SU(3) pure gauge field as given by eq. (1451) . 

Hence we find that eq. (1 127ft is the relation between the generating func¬ 
tional Z[J,7’] u ,fj u ,ri d ifidi'ils,Vs,Vh,'ijh\ in QCD and the generating functional 
Z[A, J, r) u , fj u , r) d , fjd, r} s , fj s , rjh, fjh] in the the background field method of QCD in the 
presence of SU(3) pure gauge background field A^ a (x) as given by eq. (HHft . 

Note that in QED the corresponding result is B Q 

Z[J, T ), fj] = Z[A, J, ?/, f/] (128) 


when the background field 7D(x) is the U(l) pure gauge field given by 7D(x) = <9 M u;(x). Eq. 
(135ft in QED is obtained from eq. (11281) . Note that unlike eq. (1127ft in QCD there is no 
J' in eq. (1128ft in QED because while the (quantum) gluon directly interacts with classical 
chromo-electromagnetic field the (quantum) photon does not directly interact with classical 
electromagnetic field. 

Eq. (1127ft is the main result of this paper. 

For the heavy quark Dirac field 'k(x), the non-perturbative matrix element of the type 
< 0|^(x)O n \I / (x)'l , (x / )O( l 'k(x , )|0 > in QCD is given by eq. (fT5l) if the factors O n and 0' n 
are independent of quantum fields. Similarly for the heavy quark Dirac field 'k(x), the cor¬ 
responding non-perturbative matrix element of the type < 0|'E(x)O n v k(x)'k(x / )O( l \I / (x / )|0 > 
in the background field method of QCD is given by [8] 


< 0|tf(x)O n tf(x)tf(x')O^(x')|0 > A 

= 5 O 5 5 O' 5 Z\A J 

" Mx) n Sfj h (x) S Vh (x') n 5fj h (x') [ ’ ’ 

Vsi Vh, Vh]\j=ri u =fi u =r) d =fi d =ri a =fi a =r) h =ri h =0 


Vu, Vui ’Udi Vdi Vsi 


(129) 


where the suppression of the normalization factor Z[ 0] is understood as it will cancel in the 
final result (see eq. (1133ft ). 

When the background field A im (x) is the SU(3) pure gauge as given by eq. (1451) we find 
from eqs. (fT5l) . (1129ft . (1127ft . (1124ft and (1122ft that 


< 014/ (x) On'k (x) 4/ (x 7 ) 0' n ^> (x ; ) 10 > 

=< 0|'k(x)<h(x)O n d) t (x)^(x)'k(x , )<l)(x , )O;d> t (x , )^(x')|0 


(130) 


if the factors O n and 0' n are independent of quantum fields where, see eq. (1521) . 

$(x) = exp[igT a uj a (x)\ = Ve -w f 0 °° ( 131 ) 
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Note that the creation operator aj and annihilation operator a q of the quark are related to 


the quark field via the equation 


39] 


V'(z) = \ 'TFTtMpMp ) 6 IP ' X + a\{p) v {pV P ' x ] 


spin p V V E p 


(132) 


where color indices are suppressed. Hence one finds that the quark field ip(x) or 4/ (x) 
depends on the the creation (annihilation) operator aj (a q ) of the quark but is independent 
of the creation (annihilation) operator a^ H (an) of the hadron. Similarly the gluon field 
(^“(x) is independent of the creation (annihilation) operator a^ H (a#) of the hadron. Since 
a^Of/ is independent of ^(x), 'k(x), Q 9a (x) one can perform exactly the similar steps of the 
path integral calculation as above to find from eq. (I130p that 


< 0\^(x)O n ^(x)^ H a H ^(x')O , n ^(x , )\0 > 

=< 0|#(x)$(x)O n d> t (x)d'(x)a^aHd'(x , )d>(x , )O, , l < f>t(x')'k(x / )|0 >,4 (133) 


where 4>(x) is given by eq. (11311) . 

Under non-abelian gauge transformation as given by eq. (1951) the Wilson line in QCD 
transforms as 


= u(Xf) 


Ve 


ig f dx^A“(a:)r a 


U~ L ( Xi ). 


(134) 


From eqs. (1471) and (11341) we find 
V e ~ i9 fo°° d ^- A ' a ( x + l ^ Ta = U(x)Ve~ i9 ir dXl - Aa ( x+l ^ T \ 


U(x) = exp [igT a uj a (x)} 


(135) 


which gives from eq. (11311) 

$'(x) = U(x)$(x), d>' f (x) = (x) f / —1 (x). (136) 

Hence we find that < 0|d'(x)d>(x)O n d>t(x)d'(x)a] :/ ai/d'(x , )<I>(x / )O( l < f>t(x')d'(x , )|0 in eq. 
(11331) is gauge invariant and eq. (11331) is consistent with the factorization of infrared diver¬ 
gences in QCD. 
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IX. PROOF OF FACTORIZATION IN HEAVY QUARKONIUM PRODUCTION 
IN NRQCD COLOR OCTET MECHANISM AT ALL ORDER IN COUPLING CON¬ 
STANT 


The non-perturbative matrix element < 0|T(x)O ri T(x)a^a//T(x , )O( ) T(x , )|0 > in 

QCD in eq. (I133jl is obtained from the exact generating functional in QCD as 
given by eq. m, see eq. (TT5T) . Similarly, the non-perturbative matrix element < 
0|'l r (x)O n T(x)a^a//T(x')O[ i 'h(x / )|0 >,4 in the background field method of QCD in eq. (1 1331) 
is obtained from the exact generating functional in the background field method of QCD 
as given by eq. (19UD . see eq. (I129|) . Hence we find that eq. (11331) is valid at all order in 
coupling constant in QCD. 

Note that in 4, 5] the proof of factorization is presented at NNLO in coupling constant 
and to v 2 order in relative velocity of the heavy quark-antiquark pair. This is done by 
restricting the result to v 2 order by using 


Pi = 7T +)Pr, 


n P 
2 




(137) 


where is the momentum of the heavy quark, p% is the momentum of the heavy antiquark, 
P M is the total momentum of the heavy quark-antiquark pair and p^ is the relative momen¬ 
tum of the heavy quark-antiquark pair. In the rest frame of heavy quark-antiquark pair 
p r — Mv where M is the mass of the heavy quark. 

Similarly in jfi] the proof of factorization is presented at NNLO in coupling constant 
and to all powers in relative velocity v of the heavy quark-antiquark pair. This is done by 
obtaining the result for arbitrary p± and p% without restricting to order p*. Hence in order 
to be consistent with the proof of factorization of |g| to all powers of relative velocity v it is 
necessary to present the final result for arbitrary p^ and p%■ 

It can be seen that the non-perturbative matrix clement < 
0|T(x)O n T(x)aJ / ai/'I , (x , )O( t T(x , )|0 > in QCD in eq. (j!33p is obtained from the ex¬ 
act generating functional in QCD as given by eq. (fl2|) without putting any restrictions on 
heavy quark and antiquark momenta, see eq. (1 151) . Similarly, the non-perturbative matrix 
element < 0|T(x)O n T(x)a^aH'I r (x , )O( i T(x / )|0 >a in the background field method of QCD 
in eq. (I133jl is obtained from the exact generating functional in the background field 
method of QCD as given by eq. (1901) without putting any restrictions on heavy quark and 
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antiquark momenta, see eq. (j!29j) . Hence we find that eq. (11331) is valid for any arbitrary 
momenta Pi and °f the heavy quark and antiquark respectively. This implies that eq. 
is valid to all powers in heavy quark relative velocity. 

Hence we find that eq. (I133p is valid at all order in coupling constant in QCD and to all 


Q. 


powers in the heavy quark relative velocity. 

As mentioned earlier, in NRQCD an ultraviolet cutoff A ~ M is introduced [fj. Hence 
the ultraviolet (UV) behavior of QCD and NRQCD differ. However, the infrared (IR) 
behavior of QCD and NRQCD remains same |llj. Hence the infrared behavior of the 
non-perturbative NRQCD matrix element < 0|yf K n ^(a^ H a H )^K' n x\^ > hi eq. (HJ can be 
obtained by studying the infrared behavior of the non-perturbative matrix element in QCD 
of the type < 0|T(a:)(9 n T(a:)a^a j yT(a: , )(9' 1 T(x , )|0 > where O n , 0' n are appropriate factors 
which identify the state of the heavy quark-antiquark system such as the color singlet state 
or color octet state etc.. 

We are interested in the effect of exchange of soft-gluons between the light-like Wilson 

nn 

line and the heavy quark (and/or antiquark) in NRQCD color octet mechanism |4j-|6|]. Hence 
for the color singlet S-wave non-perturbative matrix element we find from eq. (1 133ft that 

< 0|a|0 >=< x t (0)/C^(0)(aW)e t (0)ii', X (0)|0 > (138) 


at all order in coupling constant which is consistent with eq. ©■ 

When the factors O n , 0' n contain the color matrix T a we find by simplifying infinite 
numbers of non-commuting terms [see eq. ([96]) ] that eqs. (I133p and (I13ip give 

< 0\^(x)0 71ta 'i>{x)a 1 H a H ^{x , )0' na ^{x')\0 > 

=< 0\ty(x)O n ^(x)$^ ) \x)a\ { aH$ iA \x') ba ^(x')O , nia 'ir(x')\0 > A (139) 

where 

& A \x) = 'p e -i9S 0 °°™-A°{x+iX)TM° j (T^ c ) ab = -if abc . (140) 


Hence from eqs. (I139p and (1140p we find that the gauge invariant octet S-wave non- 
perturbative NRQCD matrix element which is consistent with factorization of infrared di¬ 
vergences at all order in coupling constant and to all powers in the heavy quark relative 
velocity is given by 

< 0|a|0 >=< 0|x t (0)AA^(0)# ) \0) e6 (al,aH)# ) (0) ba e t (0)A^ aX (0)|0 > (141) 
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where 


$l A \0) = Vexp[-igT {A)c f°° dM • H c (/A)], (T (A)c ) ab = -if abc . (142) 

Jo 

Note that the non-perturbative matrix element < 0\^(x)O n A(x)a^ H aH^{x')O' n ^/(x')\0 > 
in the left hand side of eq. (11331) is independent of ZO Hence all the Z M de¬ 
pendence in 4>(x) defined by eq. (1131 jl in the non-perturbative matrix element < 

in the right hand side of 

eq. (11331) is canceled by the use of background field A ,,a (x) in the expectation 
value of the non-perturbative matrix element < 0|'h(x)O rt T {x)o) H aH^{x')0' n ^> (x')|0 >a 
as dehned in eq. (11291) in the the background field method of QCD. This proves 
that the long-distance behavior of the non-perturbative NRQCD matrix element < 
0 |x t (0)ii: n ^(0)$S j4)t (0) e6 (a] ? a i /)$ z (A) (0) fea ^ t (0)^ a x(0)|0 > in eq. ( [111]) is independent of 
the light-like vector Z M at all order in coupling constant and to all powers in heavy quark 
relative velocity. 

To summarize this, we find that eq. (1 141 p . which is found by using path integral method 
of QCD, is valid at all order in coupling constant and to all powers in heavy quark rela¬ 
tive velocity. We have also shown that the long-distance behavior of the non-perturbative 
NRQCD matrix element is independent of the light-like vector Z M at all order in coupling 
constant and to all powers in heavy quark relative velocity. The eq. 03]), which is found 
by using diagrammatic method of QCD at NNLO in coupling constant and to all powers 
in heavy quark relative velocity shows that long-distance behavior of the non-perturbative 
NRQCD matrix element is independent of the light-like vector Z M at NNLO in coupling 
constant and to all powers in heavy quark relative velocity. This implies that the gauge 
invariance and the factorization at all order in coupling constant require gauge-completed 
octet S-wave non-perturbative NRQCD matrix element that was introduced previously to 
prove factorization at NNLO. 

Hence we find that eq. (]4j) is valid at all order in coupling constant and to all powers in 
the heavy quark relative velocity. 
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X. FACTORIZATION THEOREM IS A KEY INGREDIENT IN CALCULATION 
OF NRQCD HEAVY QUARKONIUM PRODUCTION CROSS SECTION 


As mentioned earlier the definition of the NRQCD heavy quarkoninm production matrix 
element from heavy quark-antiquark pair is a non-perturbative quantity which can not be 
calculated by using perturbation theory no matter how many orders of perturbation theory 
is used. From this point of view the path integral formulation (as opposed to diagram¬ 
matic methods in perturbation theory) is useful to study the properties of the NRQCD 
non-perturbative matrix element of heavy quarkonium production at all order in coupling 
constant. As mentioned earlier the only path integral formulation to study factorization 
of soft and collinear divergences at all order in coupling constant in quantum field theory 


available is by R. Tucci js] . However, the calculation of R. Tucci 


8] was exact for QED but 


was not exact for QCD. We have extended the exact path integral calculation of proof of 
factorization of R. Tucci in QED js] to proof of factorization in QCD at all order in coupling 
constant in [ 9 ] and to proof of factorization in NRQCD heavy quarkonium production at all 
order in coupling constant in the previous section. 

In this section we will show how the factorization theorem as given by eqs. (1 1301) and 
(I139p is actually a key ingredient in calculation of NRQCD heavy quarkonium production 
by using eqs. ([2]) and (|3]) where the NRQCD non-perturbative matrix element of heavy 
quarkonium production in color octet mechanism is given by eq. (1 14 1|) . 

Let us prove how the eqs. (11301) and (I139j) are key ingredients to prove eqs. (J21) and ([3]) 
to calculate the NRQCD heavy quarkonium production in color octet mechanism. Suppose 
we calculate the cross section of heavy quark-antiquark production in color octet state in 
the presence of light-like quark (or gluon). Then from eq. (I130p we find 


=< 0|t(x)$ t (x)O ri <f»(a:)T(x)'F(a; , ) < l> t (x , )O( l <f»(x / )^(a: , )|0 > . (143) 


Hence from eq. (11431) we find that the A dependence which arises due to the soft gluon 
exchanges with light-like quark (or gluon) is factorized and only appears in the gauge-links 
4>(x) in the right hand side where 4>(x) is given by eq. (I131D . Eq. (11431) implies that in the 
cross section for QQ production in color octet state at all order in coupling constant the 
infrared divergences due to the presence of light-like quark (or gluon) are factorized only to 
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the gauge links <h(x). 

Eq. (1143(1 is the exact extension of eq. (1.6) of js] of the factorization in QED. 

Hence from eq. (11431) we find that the non-perturbative matrix element of NRQCD heavy 
quarkonium production in color octet mechanism which cancels these infrared divergences 
and is consistent with the factorization theorem is obtained from eq. (I139|) and is given by eq. 
(lllip . This proves that the factorization theorem as given by eqs. (1130(1 and (1139(1 is actually 
a key ingredient to prove eqs. (J2]) and (J3]) to calculate the NRQCD heavy quarkonium 
production cross section in color octet mechanism at all order in coupling constant. 


XI. CONCLUSIONS 

Recently the proof of factorization in heavy quarkonium production in NRQCD color 
octet mechanism is given at next-to-next-to-leading order (NNLO) in coupling constant by 
using diagrammatic method of QCD. In this paper we have proved factorization in heavy 
quarkonium production in NRQCD color octet mechanism at all order in coupling constant 
by using path integral method of QCD. Our proof is valid to all powers in the heavy quark 
relative velocity. We have found that the gauge invariance and the factorization at all order 
in coupling constant require gauge-completed non-perturbative NRQCD matrix elements 
that were introduced previously to prove factorization at NNLO. 
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Appendix A: Interaction Between Non-Eikonal Current and the Gauge Field Gen¬ 
erated by Light-Like Eikonal Current in Quantum Field Theory 


From the non-eikonal part of the diagram in eq. (j7SP we find 


e 


d 4 k k v 'y I/ ^ IJt A ll (k) 
(27r) 4 2 q ■ k + ie 


= J cl 4 xJ fl (x)A fJ '(x) 


(Al) 
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where the non-eikonal current density J^(x) of the (light-like or non-light-like) charge e of 
four-momentum q M is given by 


p rOO p) 

/ d\—5 {4) (x - q\). 

2 Jo ox v 


(A2) 


Hence using eqs. 


and 


in eq. (l6Tj) we hnd that the interaction between the non- 


eikonal current and the gauge field generated by the light-like eikonal current in quantum 
field theory gives the effective (interaction) action 


s :f f [J\ = J d 4 x lv r jT dXS^(x - gA)^-^j J dX 1 6^(x - IX') 


= l 2 - 


d 4 x 


l ■ d[q ■ (x - gA 0 )] 


h 


2 f [q ■ (x - gA 0 )] 2 (/ • x) 3j 

where A 0 is the solution of the equation 


(A3) 


(x - qX 0 ) /1 (x - q A 0 ) m = 0. 


(A4) 


From eqs. (1A3h and (1A4I) we hnd that the interaction between the non-light-like non-eikonal 
current and the gauge held generated by the light-like eikonal current in quantum held 
theory gives the effective (interaction) lagrangian density 

2 e 2 (l ■ q)(q ■ x) - {I ■ x)q 2 


4 ?/(*) = ? 2 


(A5) 

(l ■ x) 3 [(g • x) 2 — q 2 x 2 } 2 

From eq. (1A5I) we hnd that the interaction between the light-like non-eikonal current and 


the gauge held generated by the light-like eikonal current in quantum held theory gives the 
effective (interaction) lagrangian density 

l\q • 0 


= 


2 (g • x) 2 {l ■ x) 3 ' 


(A6) 
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